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ABSTRACT 

We present a deep ASCA observation of a Broad Absorption Line Quasar (BALQSO) 
PG0946+301. The source was clearly detected in one of the gas imaging spectrometers, but not 
in any other detector. If BALQSOs have intrinsic X-ray spectra similar to normal radio-quiet 
quasars, our observations imply that there is Thomson thick X-ray absorption (Njj > 10 24 cm~ 2 ) 
toward PG0946+301. This is the largest column density estimated so far toward a BALQSO. 
The absorber must be at least partially ionized and may be responsible for attenuation in the 
optical and UV. If the Thomson optical depth toward BALQSOs is close to one, as inferred here, 
then spectroscopy in hard X-rays with large telescopes like XMM would be feasible. 

Subject headings: galaxies: active — quasars: absorption lines — quasars: individual (PG0946+301) — X- 
rays: galaxies 



1. Introduction 

About 10 - 15% of optically selected QSOs 
have optical/UV spectra showing deep absorption 
troughs displaced blueward from the correspond- 
ing emission lines. These broad absorption lines 
(BALs) are commonly attributed to material flow- 
ing toward the observer with velocities of up to 
- 50,000 kms" 1 . BALQSOs are probably nor- 
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mal QSOs viewed at a fortuitous orientation pass- 
ing through a BAL outflow, thus implying a BAL 
"covering factor" at least 10 - 15% in all QSOs. 
BALQSOs thus provide a unique probe of condi- 
tions near the nucleus of most QSOs. The absorb- 
ing columns typically inferred from the UV spectra 
for the BAL clouds themselves are Ah — io 20 ~ 21 
cm~ 2 (Korista et al. 1993). It has been noted, 
however, that UV studies underestimate the BAL 
column densities because of saturation (Korista et 
al. 1993, Arav 1997, Hamann 1998). BALQSOs, 
as a class, show higher optical/UV polarization 
than other radio-quiet QSOs (Schmidt & Hines 
1999, Ogle et al. 1999). Polarization studies re- 
veal multiple lines of sight through high column 
density gas (Goodrich & Miller 1995, Cohen et 
al. 1995). 

With the absorbing column densities as esti- 
mated from the earlier UV studies, we would 
have expected very little soft X-ray absorption in 
the BALQSOs. However, BALQSOs are found 
to be markedly underluminous in X-rays com- 
pared to their non-BALQSO counterparts (Breg- 
man 1984, Singh et al. 1987, Green et al. 1995). 
Green & Mathur (1996, here after GM96) studied 
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11 BALQSOs observed with ROSAT and found 
that just one was detected with a ox 15 about 2. 
BALQSOs thus have unusually weak soft X-ray 
emission, as evidenced by large a ox (>l.9. c.f. 
a ox =1.51±0.01, from Laor et al. 1997, for radio- 
quiet quasars). If BALQSOs are indeed normal 
radio-quiet QSOs, then their weak X-ray flux is 
most likely due to strong absorption. Unfortu- 
nately, due to the low observed flux, there are 
no observed X-ray spectra of BALQSOs to con- 
firm the absorption scenario, with one exception, 
the archetype BALQSO PHL5200 (Mathur, Elvis 
& Singh 1995, here after MES95). The ASCA 
spectrum of PHL5200 is best fit by a power-law 
typical for non-BALQSOs in the 2-10 keV range, 
with intrinsic absorption 2 to 3 orders of magni- 
tude higher than inferred from UV spectra alone 
(MES95). However, the PHL5200 spectrum suf- 
fers from a low signal to noise ratio, and while the 
above was a preferred fit, a model with no intrinsic 
absorption also fits the data. Recently Gallagher 
et al. (1999, hereafter G99) studied a sample of six 
new BALQSOs with ASCA, of which two were de- 
tected. G99 derived column densities of >5xl0 23 
cm~ 2 to explain the non-detections, even higher 
than the ROSAT estimates (assuming a neutral 
absorber with solar abundances unless stated oth- 
erwise) . 

How are the X-ray and UV absorbers related 
to each other? Are they both part of the same 
outflow? If so, then the kinetic energy carried out 
is a significant fraction of bolometric luminosity 
of the quasar (see Mathur, Elvis & Wilkes 1995 
for a discussion). With all QSOs likely to con- 
tain a BAL outflow, it becomes very important 
to measure the absorbing column density accu- 
rately to understand the energetics and dynamics 
of quasars. We attempt this with a deep ASCA 
observation of a typical BALQSO, PG0946+301. 

2. Observations and Data Analysis 
2.1. Observations 

We observed PG0946+301 with ASCA (Tanaka 
et al. 1994) on 1998 November 12. ASCA contains 
two sets of two detectors, SIS (Solid-state Imaging 

The slope of a hypothetical power law connecting 
2500 A and 2 kcV is defined as cx-ox = 0.384 logL pt/L 

X , SO 

that a ox is larger for objects with weaker X-ray emission 
relative to optical. 



Spectrometer) and GIS (Gas Imaging Spectrome- 
ter). The effective exposure times in SISO, SIS1, 
GIS2 and GIS3 were 72,024 seconds, 69,668 sec- 
onds, 80,910 seconds and 80,896 seconds respec- 
tively. SIS was operated in 1CCD mode with the 
target in the standard 1CCD mode position. GIS 
was operated in pulse height (PH) mode. The 
data were reduced and analyzed using FTOOLS 
and XSELECT in a standard manner (see ASCA 
Data Reduction Guide or MES95 and G99 for de- 
tails of data reduction). 

2.2. Image Analysis 

2.2.1. XSELECT Analysis 

We used XSELECT to create full and hard 
(2-9.5 keV) band images of for each of the four 
detectors. We also created combined SIS and 
GIS images. We looked for the target in these 
images displayed with SAOIMAGE. While there 
were sources seen within the GIS field of view, 
there was no obvious source seen at the target 
position in any of the four detectors. We then 
smoothed the images with a Gaussian function of 
a = 1-2 pixels. A faint source at the position of 
the target was then evident in GIS3 hard band im- 
age and a trace of a source was seen in the full GIS 
image, but not in any other image. Note that for 
a standard pointing position the target lies closest 
to the optical axis in SISO and GIS3. GIS3 is more 
sensitive in hard X-rays than SISO. The fact that 
the source is seen by eye in the GIS3 detector only 
suggests that the source is faint with flux mainly 
in the hard band. 

We extracted the total counts in a circular re- 
gion with a 3' radius centered on the source posi- 
tion. Because our source is observed to be so faint, 
background subtraction is crucial in determining 
the net source count rate, so we have done care- 
ful background subtraction using different back- 
ground estimates. Background counts were ex- 
tracted in two different ways: (1) from a source- 
free region on the detector and (2) from exactly 
the same region as the source in the blank sky 
background files provided by the ASCA guest ob- 
server facility. The significance of the source de- 
tection was therefore different for different back- 
ground estimates. For SIS, the blank sky back- 
ground is underestimated because it is available in 
the BRIGHT mode only, while the source counts 
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were extracted in the BRIGHT2 mode. So the SIS 
detections are less reliable with background (2). 
We found that the source was detected in GIS3 
and GIS3 hard band, and is marginally detected 
in SISO (2a). It was not detected in any other 
detector in either bandpass. The significance of 
detection for the source in different detectors and 
the resulting net count rate is given in Table 1. 
For non-detections, we give a 3a upper limit of the 
count rate (see G99 for exact formulation of the 
detection and corresponding count rate estimate). 

2. 2. 2. XIMA GE Analysis 

Determination of whether or not the source is 
detected is extremely important to our results. As 
an independent check, we performed image anal- 
ysis with XIMAGE (Giommi, Angellini, & White 
1997) which is designed for detailed image anal- 
ysis. The detect algorithm in XIMAGE locates 
point sources in an image by means of a sliding-cell 
method. We used detect on all of our images and 
looked for a source at the position of the target. 
Again, we found the source to be detected in GIS3 
hard band. To minimize the number of spurious 
sources detected, the threshold used by detect is 
somewhat conservative. As a result, sources with 
intensity just above the image background can be 
missed. We found that the source was detected in 
the full band GIS image if we lowered the detection 
threshold. The source was not detected in other 
detectors. These results are consistent with those 
from the XSELECT analysis discussed above. 

2.2.3. CI AO Analysis 

We applied more sophisticated wavelet-based 
techniques (Freeman et al. 2000) to provide inde- 
pendent support to the above detections. Soft- 
ware developed for Chandra Interactive Analy- 
sis of Observations (CIAO) allows us to decom- 
pose the image such that structures at different 
scales are enhanced. We analyzed the central 20' 
region of GIS3 images in both the full spectral 
range and in the harder range. Wavelet analy- 
sis of the GIS image at scales approximating the 
size of the point spread function shows that detec- 
tion of PG0946+301 is complicated by the pres- 
ence of a strong nearby source ~ 5' away. In the 
GIS hard band image, this source is significantly 
weaker, and we detect PG0946+301 at a proba- 
bility of spurious detection of 10 -4 , with a net 



count rate of (1.26 ±0.25) x 10" 3 counts s" 1 (90% 
confidence). This is consistent with the results 
discussed above. 

2.3. Column Density Constraints 

Consistency among the methods discussed 
above gives us confidence in our measurements 
and in our resulting detections in GIS3 and non- 
detections in other detectors. If the low observed 
X-ray count rate is due to intrinsic absorption, 
we can estimate the absorbing column density 
in PG0946+301. Since the source did not yield 
enough net counts in any detector to perform 
spectral analysis, we use the method discussed in 
GM96 to determine the column density. We first 
calculate the flux from the source if there was 
no intrinsic absorption. This was done using the 
observed B magnitude of the source (B = 16.0 
mag.) and assuming a ox =1.6. Redshift of the 
source (z=1.216) and the Galactic column den- 
sity (N H = 1.6 x 10 20 atoms cm -2 , Murphy et 
al. 1996) were taken into account to predict the 
2-10 keV flux in the observed band (=7.2xl0 -13 
erg s -1 cm -2 ). A power-law slope with photon 
index T = 1.7 was used. We then entered this 
model into the X-ray spectral analysis software 
XSPEC (Arnaud 1996), with normalization con- 
sistent with the expected flux and simulated spec- 
tra using SIS and GIS response matrices. The 
response of the telescope and detectors was taken 
into account as well. The column density at the 
redshift of the source was an additional parameter 
used in the simulation. If there was no intrin- 
sic absorption, then the predicted count rate was 
found to be typically an order of magnitude larger 
than the observed one. We then varied the value 
of the intrinsic absorption, keeping the normal- 
ization constant, until the predicted and observed 
column densities matched. The values of intrinsic 
column density estimated in this way are given in 
Table 2. 

This estimate of Ah depends upon T and a ox . 
Given the observed range of a ox (§l), our adopted 
value of a ox =l.6 gives conservative estimates of 
column densities. X-ray spectral slopes also vary 
among quasars. So we have estimated N# for 
r = 2.0 as well as T = 1.7. Flatter spectra result in 
even higher derived column densities. As shown in 
Table 2, even the conservative estimate results in 
Thomson thick X-ray absorption in PG0946+301, 
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i.e. Nh^IO 24 cm~ 2 . The column density esti- 
mates are consistent with the detection in GIS3 
and non detection in SISO. 

Alternatively, is it possible that PG0946+301 
(and BALQSOs in general) is intrinsically X-ray 
weak? Earlier work (GM96, G99) could not rule 
out this possibility. To test this, wc estimated the 
observed SISO hard band count rate for flux consis- 
tent with detection in GIS3 hard band, but no in- 
trinsic absorption. We find that the source would 
have been detected in SISO hard band at > 8cr 
(with T = 1.7; > 7<T with T = 2.0). So we conclude 
that the observed X-ray weakness of BALQSOs is 
due to absorption, and not due to intrinsic weak- 
ness. We cannot, however, rule out the possibility 
that the source is intrinsically X-ray weak with an 
unusual spectral shape (turning up at around 10 
keV, rest frame). It is also possible that the ob- 
served flux is only the scattered component, from 
a line of sight different from the absorbing mate- 
rial. This is unlikely in PG0946+301 which not 
strongly polarized (Schmidt & Hincs 1999). How- 
ever, if true, it again implies the existence of X-ray 
thick matter along the direct line of sight. 

3. Discussion 

We have clearly detected the quasar PG0946+301 
in our deep ASCA observation and we infer 
that there is Thomson thick X-ray absorption 
(A H > 10 24 cm~ 2 ) toward this BALQSO. The use 
of a detection, rather than upper limits, to de- 
termine the absorption is highly significant. In 
earlier work, GM96 and G99 had estimated ab- 
sorbing column densities of a few times 10 22 cm~ 2 
and 10 23 cm -2 respectively. However these were 
based on non-detections only and hence yielded 
only lower limits to the column density. A detec- 
tion provides a much stronger estimate. 

Assuming that there is indeed Thomson thick 
matter covering the X-ray source, can we infer its 
ionization state? The X-ray absorber will cover 
the optical and UV continuum sources as well, 
at least partially. If the absorber is completely 
neutral, it will result in significant HI opacity, 
which is not observed (Arav et al. 1999). If the 
absorber is completely ionized, then the opacity 
due to Thomson scattering would be the same in 
the optical, UV and X-rays (up to m e c 2 ). Thus 
this scenario by itself cannot account for the un- 



usually large values of a ox . If, on the other 
hand, the hydrogen is mostly ionized, but there 
are still some hydrogen-like and helium-like heavy 
elements, then photoelectric absorption would still 
be the dominant mechanism in X-rays. In the opti- 
cal/UV, a Thomson opacity of one would result in 
attenuation by a factor of 2.7. Such attenuation is 
inferred from polarization studies (Goodrich 1997, 
Schmidt & Hines 1999). The X-ray absorber thus 
must be at least partially ionized and may be re- 
sponsible for attenuation in the optical and UV. 

Whether the X-ray absorber has an ionization 
state overlapping the range of UV BALs and if it 
outflows with similar velocity remain outstanding 
questions. It is possible that the X-ray absorber is 
stationary, at the base of winds producing BALs. 
X-ray continuum source might be preferentially 
covered. X-ray spectroscopy is necessary to bet- 
ter probe the nuclear region in BALQSOs. For 
PG0946+301, we predict about 0.015 counts s" 1 
with the XMM PN. A reasonable spectrum may 
be obtained in about 70 ks. 
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of the US Department of Energy under contract 
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Table 1: ASCA Count Rates for PG0946+301 a 
(1CT 3 ) photons s" 1 





SISO 


SISO hard 


SIS1 


GIS2 


GIS3 


GIS3 hard 


Background* 1 


0.855(2cr) 

< 1.2 
1.29(3(7) 


< 0.74 


< 1.2 


< 0.75 


1.84(7.9(7) 


1.4(8.7(7) 


Background 2 


0.51(2a) 
< 0.74 


< 1.2 


< 0.75 


0.52(2cr) 
< 0.83 


0.48(2(7) 
< 0.63 



a. Significance of detection is given in brackets. For 
non-detections, 3a upper limit is given. For 2a 
detections, 3a upper limit is given as well. For SISO 
and GIS3, hard band count rates are given as well. 



b. With background from a source free region on the 
detector. 

c. With background from blank sky observations. 



Table 2: Column Density Constraints (10 24 atoms 
cm ) 



Detector 


r 


Detection" 


3c Lower Limit 6 


2cr Detection 


GIS3 


1.7 


0.95 


2.1 


3.3 




2.0 


0.52 


1.2 


1.95 


GIS3 hard 


1.7 


1.2 


2.55 


3.2 




2.0 


0.67 


1.55 


1.95 


SISO 


1.7 


1.4 


1.42 


1.95 




2.0 


0.9 


0.92 


1.3 


SISO hard 


1.7 




2.12 


2.73 




2.0 




1.42 


1.86 



a. If 3cr or better detection (Table 1). 

b. Upper limit on count rate gives lower limit on the 
column density. 



